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Abstract. We report in this letter the first resolved im- 
ages of the circumstellar dust disk around the Pre-Main- 
Sequence star HD 100546. These near-infrared images 
were obtained in J and short K bands with the adap- 
tive optics system ADONIS at the ESO observatory. A 
bright disk extending 2 arcsec (200 AU) far from the star 
and viewed with an inclination with respect to the line of 
sight around 50 degrees is revealed by the observations. 
Using a simplified model of light scattering, we find that 
the dust disk density peaks at a distance of 40 AU from 
the star, and the FWHM of this dense ring is typically 20 
AU. This type of disk is believed to be a denser precursor 
of (3 Pictoris-like main-sequence disks. The structure of 
the disk is compared with the best known example of the 
class of "debris disks" , the disk surrounding (3 Pictoris. 

Key words: Stars: individual: HD 100546, Methods: data 
analysis. Scattering 



1. Introduction 

The discovery of infrared excesses around an important 
fraction of main-sequence stars (Aumann et al. 1984; Plets 
& Vynckier 1999) and attributed to the presence of cool 
dust grains orbiting the star and geometrically arranged 
in a disk, has triggered a lot of studies because these disks 
may be linked to planetary formation. For a long time the 
disk around the star P Pictoris was the only example to be 
resolved both in the visible (Smith & Terrile 1984) and in 
the mid-infrared range (Lagage & Pantin 1994). Recent 
discoveries of young "debris" circumstellar disks around 
relatively "old" and isolated (i.e. not associated to any star 
forming region) stars have shown that this phenomenon 
extends towards Pre-Main-Sequence stars (cf the photo- 
metric survey by Malfait et al. 1998a). The detection and 
resolution of disk around stars such as HD 141569 (Wein- 
berger et al. 1999; Augereau et al. 1999) or HR 4796A 
thanks to high-resolution observations in the visible/near- 
infrared range (Schneider et al. 1999; Augereau et al. 1999) 
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or mid-infrared images (Koerner et al. 1998; Jayaward- 
hana et al. 1998) have shown the possibility to observe 
and study precursors to main-sequence dust disks. These 
so-called "baby-/? Pic" dust disks are the denser precur- 
sors to main-sequence debris disks. IRAS and ISO/SWS 
observations have shown that they usually produce a huge 
infrared excess (typically 250 times the infrared excess 
produced by the /3 Pic disk); some of them, as HD 100546 
for instance, showing prominent signatures of crystalline 
water ice (Waelkens et al. 1996, Malfait et al. 1998b), are 
particularly interesting targets in the visible/near- infrared 
because of high particle albedo. Observing this class of 
disks at various stages of evolution will help in finding a 
comprehensive scenario for the origin, the evolution, and 
the lifetime of the Vega phenomenon, i.e. how these disks 
form, how long they last, and how they disappear. 
We report in this paper the first images of the disk around 
HD 100546, an isolated young main-sequence star. Using 
a model of scattering by dust grains, we derive the mor- 
phology of the disk outwards 10 AU and show evidence 
for a density maximum around 40 AU. The structure of 
the disk is compared to the best known example of the 
Vega Phenomenon, the disk around /3 Pic. 

2. Observations 

We observed the HD 100546 star (mv=6.7) using the 
SHARP Il-t- camera coupled with the adaptative optics 
system ADONIS, and mounted on the ESO 3.6m telescope 
at La Silla, Chile. A pre-focal optics coronograph (Beuzit 
et al. 1997) was used to reject the direct starlight and in- 
crease the integration time in each elementary exposure. 
The mask size we used has a diameter of 0.74". J (1.25 /zm) 
and Ks (2.15 ^m) band exposures were obtained on nights 
23 and 24 of June, 1999. The seeing was quite variable the 
first night, but had decreased to a value of 0.7" at the 
time of the J and Ks observations. It was stable at a value 
around 0.85" during the second night. We spent a total 
observing time of 1550 s in the Ks band and 2600 s in the 
J band. The Point Spread Function (PSF) was frequently 
monitored thanks to interlaced observations of two ref- 
erence stars, HD 97218 and HD 101713. These reference 
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stars were used later in the reduction process to remove 
the wings of the stehar PSF and for photometric cahbra- 
tion. Observations of empty fields were also performed to 
estimate and remove the background flux which can be 
important in the Ks band. 

3. Reduction procedure 

First, standard reduction techniques were applied to the 
data : bias subtraction, flat-field correction. For each filter, 
we obtained a set of HD 100546 observations and corre- 
sponding PSFs. In spite of the use of a coronograph mask, 
the disk surface brightness is still dominated by the stel- 
lar emission at any distance from the star. In order to 
retrieve the disk image, one has to remove numerically 
the starlight wings. The rough substraction of a scaled 
PSF to HD 100546 images gives poor results because of 
slight shifts in position on the array between the refer- 
ence and the object, uncertainties on the fluxes (given 
by the literature), and a residual background (ADONIS 
bench emission, different airmasses). We developed a spe- 
cific method to get an optimum subtraction. For each cou- 
ple of HD 100546 image (Obj) and corresponding Psf, we 
have to find 3 parameters : a shift {Sx,Sy) between the 
two images, a scaling factor R, and a residual background 
Bg. These parameters are estimated by minimizing the 
following error functional : 

J = ^|Obj-5(Psf/i?,5x,(52/)-i?g| (1) 

where the S function stands for image shift. The sum is 
performed on a set of pixels (typically 1000) located in a 
region of the images where no disk emission is expected, 
and from which non accurate ones are excluded (bad pix- 
els, pixels belonging to area contaminated by diffracted 
light from the coronograph support) . The functional min- 
imum is found using a 0-order minimization algorithm 
called Powell method (Press 1996). The method was first 
checked on simulated data whose input parameters (shifts, 
scaling factor) were recovered with an accuracy better 
than 5 %. In order to evaluate the errors in the subtrac- 
tion process, test the stability of the PSF, and make sure 
that the disk obtained is not an artifact created by the 
process, the same substraction process was applied to the 
two different PSFs. Figure |l| shows the resulting images 
in the J and Ks filters of the disk and the corresponding 
residuals between the two reference stars. The position of 
the star under the mask, and therefore the disk center 
is also carefully determined in this procedure thanks to 
offset observations (out of the coronograph mask) of the 
reference stars and the object. 

4. Results 

Since the dataset in J band is more accurate (because of 
a larger integration time, better observing conditions, and 
a better sensitivity of the whole system) than in Ks band. 
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Fig. 1. Upper left : the disk seen in the J band. North 
is to the top, East to the left. The pixel scale is 0.035 
arcsec/pixel. Upper right, the corresponding residuals ob- 
tained when applying the reduction method described in 
the text to two reference stars. Lower frame : same as up- 
per frame but in the Ks filter. The "cross-like" pattern 
(with fingers orientations at 45,135,225, and 315 degrees) 
superimposed to the disk emission is produced by astigma- 
tism residuals not corrected by the adaptive optics loop, 
showing that, contrarily to J band where the correction 
residuals are still dominated by the seeing, in Ks band 
under good conditions, one can reach high Strehl ratio 
values (i.e. system residuals dominating the uncorrected 
part of the seeing) 



the geometrical parameters and the disk model are based 
on the J band image. We checked in a second step that 
the results are compatible with the Ks band data. 

4-.1- Disk geometrical parameters 

Given a disk inclination with respect to the line of sight 
which is ~ 45°, we expect to observe some deviations from 
the perfect ellipse case. Indeed, particles relatively large 
with respect to the wavelength produce non-isotropic scat- 
tering. One simple way to reduce this effect is to sym- 
metrize the disk i.e. construct an isotropic-like virtual disk 
by averaging the original disk image and its symmetrical 
image with respect to the major axis. Best ellipse fitting 
on this symmetrized images gives the following results: 

— disk position angle on the sky of 37 ± 5 degrees with 
respect to east direction. 
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— ellipse major to minor axis ratio of 1.2 implying a disk 
inclination of 50 ± 5 degrees with respect to the line of 
sight (or from the edge-on case) . 

— disk extension : « 2.0 arcsec from the star correspond- 
ing to a distance to the star around 200 AU (assuming 
a distance to the Earth of 103 parsecs (van den Ancker 
et al. 1998) and a limiting sensitivity of 5 mJy/arcsec^. 

4-2. Photometry 

Using the stars HD 101713 (B9V) and HD 97218 (G8) as 
references whose V magnitude is known and using appro- 
priate V-J and V-Ks color indices, we derived the following 
photometry for the disk in J and Ks band: a disk flux in J 
band of 0.32±0.06 Jy outside the coronograph mask and a 
maximum surface brightness of 0.5 ±0.1 Jy/arcsec^. From 
images of HD 100546 and corresponding reference star im- 
ages (HD 101713 and HD97218) obtained when shifting 
the targets outside of the mask, we can deduce a total of 
scattered flux in the J band around 1 Jy; note that we get 
also a disk-shaped emission when removing a scaled refer- 
ence to the HD 100546 image obtained outside the mask, 
but its signal to noise is very low. This result is compatible 
with the infrared excess computed considering a Kurucz 
model of the star (B9V spectral type). Consequently, the 
disk flux outside the coronograph mask represents roughly 
30% of the total disk emission. 

In the Ks band, we find an emission of 0.25±0.1 Jy outside 
the coronograph mask and a maximum surface brightness 
of 0.25 ±0.1 Jy/arcsec2 

4-3. Model and observation fit 

In order to derive physical parameters of the disk, we built 
a numerical model of starlight scattering by an optically 
thin disk made of dust particles. Assuming that no mul- 
tiple scattering is occurring in the disk, and since the ex- 
act physical properties of the particles are poorly known 
(composition, shape, size distribution), we chose to use a 
global, single-parameter, scattering phase function. The 
most practical one, although not realistic, is the Henyey- 
Greenstein phase function (Henyey & Greenstein 1941). 
The disk is assumed to be seen with an inclination de- 
rived above, and to have a proper thickness similar to the 
one of P Pic parameterized following Artymowicz et al. 
1989. The disk midplane density is described by a scat- 
tering area which follows a series of broken power laws. 
We start with a first guess deduced from the radial sur- 
face brightness, and we iterate on the parameters until we 
obtain a satisfactory fit along the four directions defined 
on figure |[ 

4- 4- Fit results 

The best fit model provides the following parameters (one 
must keep in mind that they are model-dependent) : 
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Fig. 3. The normal optical thickness of the disk deduced 
from model fitting of the data (diamonds). We have over- 
plotted (stars) the normal optical thickness found in the 
case of the disk around (3 Pictoris (Pantin, Lagage & Arty- 
mowicz 1998) 

— The Henyey- Greenstein phase function parameter 
found is g = 0.2 corresponding to dominant spherical 
particles with size around 0.1 micron. 

— A radial normal optical thickness fitted by a series of 
broken power laws whose exponents are : 0.6 in the 
range [10,45]AU, -0.4 in the range [45,70]AU, -1.35 
in the range [70,100]AU, -5 outwards from 100 AU, 
see figure ^ 

— An exponentially decreasing vertical profile, and an 
opening angle around 0.1 radian. 

— A total scattering area of 10^^ cm^ and a dust mass of 
0.02 Me. 

— An approximate disk surface density at maximum 
around 3x 10^* particles/m^ assuming single-sized par- 
ticles of 0.1 /im, or an equivalent maximum scattering 
area of the order of the unity/m^. 

5. Discussion 

The radial surface brightness distribution shows that the 
disk density peaks at a distance of 40 AU from the star cor- 
responding to the edge of the planetary zone in our solar 
system. The normal optical thickness found is more than 
100 times larger than in the case of the /? Pic disk. Simul- 
taneously, the presence of sub-micronic particles sizes, well 
below the radiation pressure cut-off in sizes, suggests that 
huge quantities of very small particles are produced con- 
tinuously by evaporation or collisions of planetesimals. We 
are probably witnessing a analogue of the young Kuiper 
belt with planetesimals undergoing violent agitation. If 
the dust particles are the result of collisions, a high rate is 
required to produce the observed disk. This implies that 
gravitational perturbations occur frequently in this disk, 
maybe because of giant planets. If the disk is replenished 
by planetesimal evaporation as suggested by Lecavelier des 
Etangs (1998) in the case of the (3 Pic disk, 10^ planetesi- 
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Fig. 2. Fit of the data profiles. Up- 
per plot, full lines : the modeled pro- 
files, stars and diamonds represent 
the data profiles 1 and 3 defined 
on the disk image inserted. Lower 
plot, dashed line : model of pro- 
file 4 (dominant forward scattering), 
plain line : model of profile 2. As 
seen on the lower plot, the scattering 
asymmetry is well reproduced by an 
anisotropic phase function (Henyey- 
Greenstein, see the text) with g — 
0.2, corresponding to a particle size 
around 0.1 /xm. 
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mals (40 km size) whose orbits were recently perturbed by 
planet migration are required to maintain the disk. Some 
of them which orbits are gravitationally perturbed might 
be the origin of the casual redshifted absorption profiles 
observed in some of the spectra of HD 100546 (Grady et 
al. 1997), a phenomenon highly similar to what is observed 
in the system of (3 Pictoris (Lagrange-Henry et al. 1987, 
Beust et al. 1996). 

6. Conclusion 

Thanks to coronograph adaptive optics observations we 
could detect the disk surrounding the PMS star HD 
100546. This detection of a disk which is believed to rep- 
resent a class of precursors to (3 Pic-like disks adds a new 
clue in the understanding of the formation and evolution 
of these debris disks. Other imaging observations of disks 
at various stage of evolution are necessary to sketch more 
closely the possible scenario of such a disk creation and 
evolution. This should be achieved by next generation in- 
struments such as the VLT/NAOS (adaptive optics) and 
VLT/VISIR (mid-infrared) instruments or the NGST. 
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